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Abstract. In this note, we discuss the problem of quark-antiquark pair production in the framework of the 
color glass condensate. The cross-section can be calculated in closed form for the case of proton-nucleus 
collisions, where the proton can be considered to be a dilute object. We find that fc^-factorization is broken 
by rescattering effects. 

PACS. 11. 80. La Multiple scattering - 11.15.Kc Classical and semi-classical techniques 



1 Introduction 

The Color Glass Condensate (CGC) [1-3] provides a fra- 
mework for studying hadronic collisions at high energy. 
Following the ideas developed in the McLerran- Venugopa- 
lan Model [4-7] , it divides the partonic degrees of freedom 
into hard color sources p - partons with a large momentum 
fraction x that travel at almost the speed of light and that 
can be considered as frozen over the typical interaction 
times - and classical color fields . The latter represent 
the small x modes, which have a large occupation number 
and thus can be described classically. The two types of 
modes are related by the fact that the classical fields obey 
the Yang-Mills equations, with a source term given by the 
hard color sources. 

The hard sources are random variables (the partons 
at large x come in a different configuration in each colli- 
sion), whose statistical distribution is described by a func- 
tional density W Xo [p] . The subscript x in this functional 
is the separation between the degrees of freedom that 
arc described as sources and those that are described as 
fields. Upon changing this boundary, the functional W Xo 
must change according to a renormalization group equa- 
tion, known as the JIMWLK equation [8-19]. The initial 
condition for this evolution equation is in principle non- 
perturbative. However, in the case of a large projectile 
like a nucleus, it has been argued by McLerran and Venu- 
gopalan that a good model for this initial condition is 
given by a Gaussian [4-7] . Such a Gaussian should remain 
a reasonable description of the hard sources in a nucleus 
as long as Xq does not become too small. 

In a collision between two hadronic projectiles, each 
projectile is described by its own color source pi.2, which 
correspond to partons moving at the speed of light in op- 
posite directions. The first step in studying this collision 
process is to solve the Yang-Mills equations in the presence 



of these two source terms. So far, this has been achieved 
analytically only in the case where at least one of the 
sources is assumed to be a weak source and treated per- 
turbatively (i.e. the solution of Yang-Mills equations is ex- 
panded order by order in powers of this weak source) [20, 
21]. This situation of course prevails in the collision be- 
tween two small, non saturated, projectiles, like two pro- 
tons at not too small x. It is also relevant in the case of a 
small and a large projectile, like proton- nucleus collisions, 
or deuteron- nucleus collisions as performed at RHIC. The 
solution of Yang-Mills equations in this asymmetrical sit- 
uation will be reviewed in section 2. The case of nucleus- 
nucleus collisions, where none of the color sources can be 
considered to be weak, is at the moment out of reach of 
analytic calculations, but can be solved numerically [22- 
26]. 

Once the classical gauge fields are known, physical 
cross-sections can be calculated. It is straightforward to 
obtain the gluon multiplicity 1 at the classical level, since 
it is simply a matter of doing a Fourier decomposition 
of the gauge fields. In order to obtain the multiplicity of 
produced quark-antiquark pairs, one needs the retarded 
propagator of a quark in the presence of the previously de- 
termined classical gauge field. Again, in situations where 



In the case where both sources are strong, and thus the 
classical gauge fields are also strong fields, it is much sim- 
pler to evaluate inclusive quantities like multiplicities instead 
of cross-sections for more exclusive processes. Indeed, the latter 
involves the calculation of time-ordered amplitudes in the pres- 
ence of the classical gauge field. The fact that this gauge field is 
time dependent in a collision process renders this problem very 
difficult in general due to the presence of vacuum-to-vacuum 
diagrams [27]. On the contrary, the calculation of multiplici- 
ties only involves retarded amplitudes in the classical field, for 
which there are no vacuum-to-vacuum diagrams. 
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the background gauge field is only known numerically, the 
quark propagator can only be computed numerically as 
well [28-30]. For collisions between a small and a large 
projectile however, one can obtain the quark production 
cross-section analytically, as was already the case for the 
classical color fields. The result of this calculation will be 
presented in section 3. 

Finally, as a last step involved in calculating an observ- 
able relevant for hadronic collisions, one must perform the 
average over the configurations of the hard color sources, 
weighted by a factor W[p] for each projectile. This aver- 
age cannot be performed analytically in general, except in 
some special cases like when W[p] is a Gaussian functional. 
The source averages that we need in the expression of the 
quark production cross-section will be given in section 4. 

At this point, the quark production cross-section is 
expressed as a multi-dimensional integral that cannot be 
further simplified analytically. One observes that /^-facto- 
rization is broken by re-scattering effects, contrary to what 
happened in the simpler case of gluon production. Some 
implications of this formula are discussed for single quark 
production and for quarkonium production in section 5. 



2 Gauge fields 

When studying the collision of two hadrons in the frame- 
work of the Color Glass Condensate, the first step is to 
determine the classical color field created in this collision 2 . 
One has to solve the Yang-Mills equations, 



[D li ,F^\ = r 



(1) 



For these equations to be compatible with the constraints 
provided by the Jacobi identity, the color current J u that 
appears in the r.h.s. must be covariantly conserved: 



[D v ,r\ = 0. 



(2) 



This condition implies that the current in general receives 
corrections to all orders in the color sources. At lowest 
order in the sources, it is given by 



is to obtain the solution to first order in p\ (the source 
describing the proton) and all orders in p 2 (the source 
describing the nucleus). In the Lorenz gauge, this solution 
has a nice and compact form, which reads, in Fourier space 
[21], 



' q z + iq+e 



+C^(q) [V(k 2± ) (2nf5(k 2± )] .(.I) 



k 2 



In this formula, the first term is the color field of the 
proton alone: 



A£(q) = 2ng5» + 5( q -)?^ 



(6) 



In the second term, fcu_ is the momentum coming from 
the proton and k 2 ±, defined as k 2 ± = q ± — ki±, is the mo- 
mentum coming from the nucleus. The 4-vectors 3 and 
C£ have been introduced in [21], and U, V are the Fourier 
transforms of Wilson lines in the adjoint representation of 
SU(N): 



U (xj_) = V+ exp 
V{x±) = V+ exp 



ig f dz+A-(z+,x ± )-T 
i\ J + °° dz+A-(z+, X± )-T 



(7) 



where A~ is the gauge field of the nucleus alone. The T a 
are the generators of the adjoint representation of SU(N) 
and V+ denotes a "time ordering" along the z + axis. The 
peculiarity of this result is that it contains a Wilson line 
(V) with an unusual factor 1/2 in the exponent. Such a 
Wilson line must be an artifact of the gauge we use, and 
its cancellation from physical quantities is a non-trivial 
test of the gauge invariance of the final result 4 . 



J v {x) = gS v+ 5(x-)p 1 (x ± ) + g5 v -5{x + ) P2 { Xl _) , (3) 

where p\ and p 2 are the distributions of color charge in 
the hadron moving to the right and to the left respectively. 
The previous three equations must be supplemented by a 
gauge condition. In our work, we have used the Lorenz- 
covariant gauge 

= . (4) 

Thus far, the solution of this set of equations is not 
known to all orders in both p\ and p 2 (although numerical 
solutions have been obtained). The best one can achieve 



2 An exception to this is the case where one treats one of the 
two hadrons using the standard collinear factorization instead 
of a description based on the CGC [31-36]. In this case, one 
needs only the gauge field produced by the other hadron alone 
(which is of course much simpler). 



3 Quark production 

3.1 Pair production amplitude 

Once the gauge fields have been obtained, the probabil- 
ity of producing a quark-antiquark pair in the collision is 



3 They are related to Lipatov's effective vertex by + 

4 The gauge field produced in pA collisions has also been 
determined in the Fock-Schwinger gauge x + A~ + x~A + = 
[20] and in the gauge A + — [37]. It has been checked explic- 
itly that although the gauge fields are not the same in these 
other gauges, they all lead to the very same expression for the 
multiplicity of produced gluons. 
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given by 5 : 



d 3 q 



16tt 3 



d 3 p 



(27r) 3 2w F 



\u(q)T(q,-p)v(p)\ 2 , (8) 



where u q = y/ q 2 + m 2 is the on-shell energy of a quark, 
and T(q, —p) the time-ordered propagator (amputated of 
its external legs) of a quark in the presence of the previ- 
ously obtained color field, with incoming momentum — p 
and outgoing momentum q. This probability must be inte- 
grated over the impact parameter of the collision in order 
to obtain the pair production cross-section. 

As explained in [38], there is one technical complica- 
tion when using the gauge field given in eq. (5) in order to 
calculate the pair production amplitude. Indeed, it turns 
out that one of the terms in the gauge field is proportional 
to S(x + ) in coordinate space, which means that such a 
term allows the vertex where the pair is being produced 
to be located inside the nucleus. Practically, this means 
that this term must be considered separately, because its 
contribution can only be calculated by smearing out the 
thickness of the nucleus. When this is done properly, one 
obtains a pair production amplitude in which the spurious 
Wilson line V does not appear any longer: 

... \ 2 fd 2 k 1± d 2 k± p P:a (k 1± ) f 2 2 

MAq,P) = 9 j^rj^Y JL 1 ^jd *xd y x 

xe ik ± -x ±e i(p ± +q ± -k ± -k 1± )-y ± 

Xu(q){T q g(k lx ,ks_)[U(xs_)t a W(y x )] 
,(k 1± )[t b U ba (x x )]}v(p) , 



(9) 



where we denote: 
T qq (k 1± ,k ± ) = 



If: + m)7 (4 — ft — fti + m)^ 



2p+[{q ± -k x ) 2 +m 2 } +2q+[{q ± ~k x -k 1 _ 



T g (k 1± ) EE 



(P + qf 



(10) 



In eq. (9), U is the same Wilson line as U, except that 
it is now in the fundamental representation of the gauge 
group. 

3.2 Pair cross-section 



both projectiles, and integrating over the impact parame- 
ter of the collision in order to obtain the cross-section for 
pair production. This gives 



d 2 q ± dy q d 2 pj_dy p 
a 2 N 



8it 4 (N 2 - 1) 



S{p x + q x -k 1A _-k 21 _) 
h 2 h 2 



< I tr 

/ k i ,k' , 



ki±.k 2 ± 

{+m)T qq -y-m)T;'] cj ) f^(k 2± -k ± ,k' ± ) 



tr 



+ m)T q S-m)T: ^f' ff (fe 2 x; k±) + h.c. 



-m)T g (p'-m)T* 



0r( fe 2±)}^ P (fci±) , (ii) 

The function (p p , and the various A 's are correlators of 
color sources, defined as follows [38]: 

<f>™(lx) = ^-S^ / e a ^ (U(Q)Ui(x ± )) aa , 



g 2 N 



4>f' 9 {lv,kA_) 



2n 2 R 2 A l\ 
g 2 N 



e 



i(k± ■x± + (l±-k±)-y ± ) 



bf^(l ± ;k ± ,k' ± ) = 



xtr{U(x ± )t a U^y ± )t b U ba (0)j , 
2nl 2 



g 2 N 



e i(k±-x±-k' ± -u±) e i(l±-k±)-y ± e -i(l±-k' ± )- 
x±,y ± ,u±,v± 



xtr (U(x ± )t a U\y ± )U(v ± )t a U^(i 



(12) 



The momentum assignments for the various distributions 
we have defined are illustrated below: 



<P>*(1 ± ) EE 



bf 9 (l±;kj 



ill t\i 



At this point, it is a simple matter of squaring this am- 
plitude, performing the average over the color sources in 



5 As shown in [27], this formula is in general incomplete 
because it does not include the contribution of vacuum-to- 
vacuum diagrams, which is in principle required by unitarity. 
However, in our approximation where only the leading order 
in the source pi is kept, this correction can be safely neglected. 
Note also that in this approximation, it would be equivalent to 
use the retarded quark propagator, instead of the time-ordered 
propagator. 



/ -h tf L* rZ tlf 

'A "lis I i I Bl t ± K X 



gl^ . (13) 



W(lr,kxM= W ^ 

A ( ir ? i f if } 

It is not possible to write this cross-section in terms of a 
single distribution describing the nucleus, and therefore a 
strict fcj^-factorization is impossible to achieve 6 . However, 



6 Let us recall that in the case of gluon production, the cross- 
section can be expressed in terms of the function (j>^ 3 only [40, 
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a relaxed form of fc^-factorization remains true, termed 
"non- linear /c_i_-factorization" in [41,42], provided one in- 
troduces three distinct functions describing the content 
of nucleus. The three </> A 's are not entirely independent, 
since they are related by some sum rules [38], but they 
nevertheless probe different correlations among the color 
sources contained in the nucleus. 



3.3 Single quark cross-section 

It is also possible to integrate out the phase-space of the 
antiquark in order to obtain the inclusive single quark 
cross-section 7 . One obtains [38] 



are of interest in assessing issues such as the magnitude 
of the breaking of /cj_-factorization, or in order to make 
phcnomcnological predictions in this kinematical domain. 
The expressions of the various correlators have been de- 
rived in [38]. Let us denote 

r(x± - y ± ) = g 4 / [G (x± - z_l) - G (y ± - z ± )} 2 , 

(17) 

where Go is the free 2-dimensional propagator 



G (x A _ - z ± ) 



d 2 k± e ik±-( X ±-z ± ) 



(18) 



d 2 q ± dy q 8ir 4 {N 2 
x jtr 



i) J p + J kl± , 



i 



+m)T q $-m)T^\ -f<t>™{k 2± ) 



(2^) 2 k\ 
The simplest ones are the 2-point correlators 

(U(0)Ui(x ± )) aa = (N*-l)e-^ r ^ , 

tr (u(0)U*(x±)) = Ne-^^A r ^ x ^ . (19) 



Jk ± 



+ I iv 

+tr 



(4+m)T qq -W-m)T; 



»(fc 2± ;fc ± )+h.c. 



(4+m)T g ^-m)T;]r j : 9 (k2±)}^ P (k 1± ) , (14) 

A new distribution appears in the single quark cross-sec- 
tion, defined as: 

*T C-0 = f x± ^' x ^ (u(o)uH^)) ■ (i5) tr ( 



Since the 3- and 4-point correlators have quite complicated 
expressions, we will only quote here the simplified results 
for the large N limit 8 



n 2 .[r( x± )+r( Vl _)] 



Once again, strict fc^-factorization is broken: one needs 
three different correlation functions in order to write this 
cross-section. 



4 Correlators of Wilson lines 

It is in general extremely difficult to calculate the correla- 
tors of Wilson lines that appear in the previous formulae 
for the pair and single quark cross-sections. In principle, 
one would have to solve the JIMWLK equation, which is 
a question that so far has not received a full (numerical 
or otherwise) answer, despite some encouraging attempts 
[43]. 

One can however calculate these correlators analyti- 
cally in the McLerran-Venugopalan model where the dis- 
tribution of color sources in the nucleus has the Gaussian 
form 

f Pa(x ± )p a (x_ 
J x ± 



tr(u(x ± )t a U t (y ± )t b U ba (0) i 
= NC f e 

= NC f e~ I * r ' J ' 2 A lr( - x± ~ u± ' >+r ( y± ~'" ± ^ . (20) 



W[p] = exp 



(16) 



Since this model is quite relevant at moderate values of x 
for large nuclei, the closed expressions one obtains from it 



One sees that the 3- and 4-point correlators are expressed 
entirely in terms of two-point correlators in the large N 
limit. This simplifies the computation enormously. 



5 Discussion 

5.1 Generalities on fc^-factorization 

The two terms in eq. (9) have a simple interpretation. The 
first term, containing two Wilson lines in the fundamen- 
tal representation, corresponds to a process in which the 
quark-antiquark pair is produced before the collision with 
the nucleus. The pair then goes through the nucleus and 
scatters off the color charges present in the nucleus. The 
second term, which contains only one Wilson line in the 
adjoint representation, corresponds to a process in which 
the gluon coming from the proton goes through the nu- 
cleus and produces the quark-antiquark pair only after the 
collision with the nucleus. 
Thanks to the identity 



20,21]. Therefore, even though this function is not the canon- 
ical unintegrated gluon distribution of the nucleus, one has a 
cross-section which is exactly fcx-factorizable. 

7 Single quark production has also been studied by Tuchin 
for Gaussian correlators. Unfortunately, his results are in co- 
ordinate space and cannot be easily compared to [39] 



u(x ± )ru*(x±) = t b u ba ( X j 



(21) 



It is possible, albeit quite cumbersome, to evaluate the 3- 
point function numerically without doing the large N approx- 
imation [44]. 
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between Wilson lines in the fundamental and adjoint rep- 
resentations, one sees readily that ^-factorization is re- 
covered in the limit where the difference x± — y ± between 
the transverse coordinates of the quark and the antiquark 
is neglected. In other words, fcj_-factorization should be 
a valid approximation in any kinematical regime where 
the typical transverse size of the pair is small. This is ex- 
pected to be the case if the quark mass is large, or if the 
transverse momenta in the final state are large. 

There is another way to reach the same conclusion 
from eq. (9). If the dependence of the factor T qq -(ki±, k±) 
on k±_ is neglected, then the integration over fc^ in the 
amplitude produces a S(x± — y_|_). This automatically en- 
sures /c^-factorization. Since the typical k± is of order Q s , 
it is possible to neglect k± in the factor T qq if there is some 
other momentum scale (either m or q ± ) much harder than 

Qs- 

5.2 Quark production 

The simplest physical observable to look at is the cross- 
section for single quark production, obtained by integrat- 
ing out the phase-space of the anti-quark. This partonic 
cross-section can be converted into the cross-section for D 
or B mesons production, by convolution with the appro- 
priate fragmentation function. Since quite a lot of phe- 
nomenology has been done in the framework of fcj_-fac- 
torization, it is particularly interesting to investigate the 
importance of the breaking of this factorization, and its 
dependence on parameters such as the quark mass, the 
quark transverse momentum, and the saturation scale in 
the nucleus. 

A numerical computation of eq. (14) is under way, and 
detailed results will be reported elsewhere [44]. One can 
nevertheless point to the following qualitative trends: 

— The magnitude of the breaking of /c_L-factorization de- 
creases as the quark mass increases. Indeed, since the 
terms that break ^-factorization correspond to extra 
rescatterings, it is natural that massive quarks are less 
sensitive to these effects than light quarks. 

— The magnitude of the breaking of fc^-factorization is 
maximal for a transverse momentum q± ~ Q s of the 
quark, where Q s is the saturation scale in the nucleus. 
One recovers ^^-factorization when the quark trans- 
verse momentum becomes much larger than all the 
other scales. 

— If Q s remains smaller or comparable to the quark mass 
and transverse momentum, the corrections due to the 
breaking of /^-factorization enhance the cross-section. 
This is interpreted as a threshold effect: having more 
rescatterings tend to push a few more QQ pairs just 
above the kinematical production threshold. 

— If Q s is large compared to the mass of the quark, then 
the corrections due to the breaking of k± -factorization 
tend to reduce the cross-section at small transverse 
momentum. Since the typical momentum transfer in a 
scattering is of the order of Q s , it is indeed more diffi- 
cult to produce light quarks with a small momentum 
if they scatter more. 



In addition, it will be interesting to investigate the effect 
of this breaking of /c^-factorization on the Cronin effect 
and its rapidity dependence, in order to see whether the 
Cronin effect for heavy quarks follows the same pattern 
as for gluon production [21,45,46]. 

5.3 Quarkonium production 

Quarkonium production in high-energy heavy ion colli- 
sions is a key observable for diagnosing Quark Gluon Plas- 
ma formation [4 7]. The anomalous suppression of the J/tp, 
compared with the baseline of the nuclear suppression, 
was indeed observed at CERN-SPS[48]. In the experimen- 
tal analysis, the normal dissociation of the J/ip state by 
the nucleons is assumed to be independent. At higher ener- 
gies, however, the nuclei in their center-of-mass frame are 
Lorcntz-contracted to be narrower than the characteristic 
size of the quarkonium. Thus one expects the coherence 
between the interactions with the nucleons to be more im- 
portant. In the rest frame of a nucleus, on the other hand, 
the internal motion of the heavy quarks is frozen during 
the interactions. 

A model of the quarkonium projected on a target nu- 
cleus at high energy is considered in [49,50], which is for- 
mulated as an eikonal propagation of the bound heavy 
quarks in the random color background fields in the nu- 
cleus. The multiple interactions with the color gauge fields 
result in the random walk of the quarkonium state in the 
momentum and color spaces. The survival probability of 
the bound state, the overlap of the initial and final pair 
state, is shown to be damped like a power-law rather than 
exponentially in the effective target thickness L. This phe- 
nomenon is dubbed "super-penetration" in the QED case. 

In order to understand the baseline of the normal nu- 
clear suppression in the high-energy heavy ion collisions, 
it would be crucial to investigate the nuclear effect on 
the production process of quarkonium rather than the at- 
tenuation of the formed resonance [5 1,52]. In this context, 
the analytic expression for the quark production ampli- 
tude in the proton-nucleus collision within the MV model 
will provide a good starting point for the further study of 
quarkonium formation. The multiple scattering effect in 
the nucleus target is appropriately included in the frame- 
work, and it will suppress the probability for producing the 
quark-pair with small relative momentum. We speculate 
that the formation probability of the small-momentum 
pair will be suppressed in some power of the target size, 
and not exponentially, reflecting the random walk nature 
of the multiple scatterings. 

6 Conclusions 

The calculation of the production quark-antiquark pairs in 
the framework of the Color Glass Condensate, in the case 
of proton- nucleus collisions, shows unambiguously that 
fc^-factorization is broken. This is to be contrasted to 
what was observed for gluon production where this fac- 
torization is preserved, despite the presence of important 
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rescattering effects in the nucleus. In the case of quark pro- 
duction, one needs more correlators in order to describe 
the nucleus. In particular, the cross-section depends on 
correlators of three and four Wilson lines. This opens the 
possibility of obtaining more detailed information about 
correlations among color charges, as a function of rapid- 
ity. Conversely, experiments can test in detail CGC pre- 
dictions for these higher point correlators. 
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